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ABSTRACT 

In this photometric study of the stellar association LH 95 in the Large Magellanic Cloud (LMC) we 
focus on the pre-main Sequence (PMS) population in order to construct, for the first time, the sub- 
solar initial mass function (IMF) in the LMC. The basis for this investigation consists of the deepest 
photometry ever performed in the LMC with the Advanced Camera for Surveys (ACS) on-board the 
Hubble Space Telescope {HST). We improve our catalog of point sources obtained with ACS (Paper I). 
We carry out a Monte Carlo technique to subtract the contribution of the general field of LMC. We 
isolate the central region in the observed area of the association, as the part characterized by the 
highest concentration of PMS stars. We analyze the reddening distribution of the system in order to 
obtain the mean value of extinction and we study the mass function of its field-subtracted population. 
For this purpose, we introduce a new set of evolutionary models, derived from the calculations on 
the evolution of PMS stars by Siess and collaborators. We use these models with our observations of 
LH 95 to derive the IMF of the system. This mass function is reliably constructed for stars with masses 
down to ~ 0.43 Mq, the lowest mass ever observed within reasonable completeness in the Magellanic 
Clouds. Consequently, its construction offers an outstanding improvement in our understanding of 
the low-mass star formation in the LMC. The system IMF of LH 95 shows a definite change in its 
slope for masses M <1 Mq, where it becomes more shallow. In general, the shape of this IMF agrees 
very well with a multiple power-law, as the typical Galactic IMF, down to the sub-solar regime. The 
change in the slope ("the knee") of our IMF at ~ 1 Mq also agrees with the average Galactic IMF. 
As far as the slope of this system IMF is concerned, it is found to be somewhat more shallow than the 
corresponding "classical" Galactic IMF in the sub-solar regime, probably due to unresolved binarity, 
while for stars with M <; IMq it becomes slightly steeper. We do not find significant differences in 
the shape of the overall IMF of LH 95 from that of each of the three individual sub-clusters of the 
association. This clearly suggests that the IMF of LH 95 is not subject to local variability. 
Subject headings: stars: formation — stars: pre-main-sequence — open clusters and associations: 
individual(LH95) — Magellanic Clouds — initial mass function 



1. INTRODUCTION 

The Magellanic Clouds, the closest undisrupted dwarf 
galaxies to our own galaxy, host numerous young stel- 
lar associations. The study of the low-mass populations 
of these systems provide important improvements in un- 
derstanding extragalactic star formation. The investiga- 
tion of the star formation processes, and the Initial Mass 
Function (IMF), together with the time dependency of 
the formation events, are key points in the characteriza- 
tion of stellar populations in any kind of concentrations, 
from small clusters to entire galaxies. 

Regarding the IMF, several studi es have been carried 
out in the past for the G alaxy (e.g. Salpeteilll955t IScalol 
119861 lKroupal[200ll : [Chabrier, ,2003). describing the ini- 
tial numbers of massive stars (M > 8 Mq), stars of 
intermediate mass (1 Mq < M < 8 Mq), low-mass 
stars (0.08 Mq < M < 1 M©) and brown dwarfs 
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{M < 0.08 Mq). The findings of such investigations 
support the hypothesis of the universality of the IMF, 
in the sense that the average IMF measured in a diver- 
sity of galactic syster ns does not show noticeable vari- 
ations (jKroupal l2002l ). This statement turns out to 
be rather certain for the intermediate- and high-mass 
stars, while for low-mass stars there could be a depen- 
dence of the IMF slope on metallicity, in the sense that 
metal-rich environments tend to produce more lo\y-mass 
stars t han metal-poor systems (jPiotto fc Zoccalil 119991 : 
I Re vie fc Robin 2001j), although several IMF measure- 
ments are needed to confirm such a systematic effect. In 
low-mass studies there are several uncertainties present. 
For example the stud y of the IMF in metal-p oor galac- 
tic globular clusters (jPiotto fc Zoccalil Il999( ) is signifi- 
cantly limited by the corrections for dynamical evolution 
effects. Specifically, phenomena such as mass segrega- 
tion and unknown binarity fraction change dramatically 
the presently measured mass distribution with respect 
to the one at the time of the formation of these systems 
(iMarks et a l. 2008). 

Under these circumstances, ideal environments for the 
search of a metallicity dependence of the IMF are the 
Magellanic Clouds. Their stellar populations are charac- 
terized by metallicities of [Fe/H] ~ 2.5 to 5 times lower 
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than in the disk of the Milky Way (iLuck et al.lll998D. but 
have a similar star formation rate (SFR) |Westerlundl 
|1997() . Furthermore, the low metal abundances of the 
Magellanic Clouds suggest that their environments — as 
well as in other galaxies of the Local Group — are closer 
to the ones at the early ages of the uni verse, when th e 
peak of star formation occurred {z ~ l.S. lPei et al.lll999l ). 
Moreove r, the lower dust-to-gas ratio in the M agellanic 
Clouds (iKoornneej 119821 : iBouchet et al.l 119850 . in com- 
parison with the Milky Way, assures a lower extinction, 
reducing the bias introduced by differential reddening, a 
well known limitation in the study of young stellar pop- 
ulations. 

Most of the star formation occurs in dense cores 
of Giant Molecular Clouds, where one or more 
large stellar concentra tions, known as OB associations 
(jAmbartsumian 1 11947[ ). are formed. However, most of 
the star formation studies in the Milky Way are car- 
ried out in smaller fields, like the young Taurus, Lupus 
and Chamaeleon star- forming regions. Although it has 
been shown that galactic OB asso ciations host popula- 
tions of faint, low-mass PMS stars (iPreibisch et alll2002l : 
iSherrv et 311120041 : iBriceno etld] l2007f ) their analysis is 
strongly limited by the contamination of background 
and foreground evolved populations, requiring detailed 
measurements of proper motions, time-consuming spec- 
troscopy and multi-epoch photometry to define their 
membership. On the other hand, the study of OB as- 
sociations in the LMC is less subject to contamination 
by field po pulations, due to the limited distance spread 
(|Caldwell fc Coulson 1986; Cole 1998). However, in or- 
der to identify the faint red PMS stellar component of 
LMC associations, photometry with instruments of high 
sensitivity and angular resolution is required, due to the 
larger distance from us. 

Recent findings from the Hubble Space Telescope 
(HST) confirm the presence of PMS stars in OB as- 
sociatio ns of the Magellanic Clouds. iGouliermis et al.l 
()2006aD studied this case in the LMC association LH 52, 
using WFPC2 observations in the V- and /-equivalent 
bands, and they discovered a population of objects in 
the color-magnitude diagram (CMD) consistent with 
T-Tauri stars. Subsequently, the Advanced Camera 
for Surveys (ACS) has been used to perform similar 
studies in other associations of the Magellanic Clouds, 
such a s NG C 346 (Nota ct al. 2006; Gouliermi s et all 
I2006bt Hennckemp er et all 120081 : ISabbi et al.ll200l and 
NGC 602 (Schmalzl et al.l 12001 ) in the SMC. However, 
these studies could not characterize the sub-solar mass 
function, because of the insufficient detection of such 
stars. 

In the first part of our study ACS observations en- 
abled us to discover a rich population o f PMS stars in the 
star-f orming region LH 95 in the LMC (| Gouliermis et all 
120071 from hereafter Paper I). In this paper we focus 
on the IMF of these stars, which is reliable for masses 
down to about 0.3 Mq, considering that our photometry 
allowed to perform precise completeness corrections for 
stars with masses down to this limit. In §[5] we describe 
our photometry and the completeness of our data. In 
§ [3] we apply the decontamination of the observed stel- 
lar sample in the area of the association LH 95 from the 
contribution of the general LMC field for the identifi- 
cation of the stellar members of the system. We also 
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Fig. 1.— F555W - F8UW, F8UW Color-Magnitude Diagram 
of the entire LH 95 region, according to the updated photometric 
catalog, which includes about 900 stars more, recovered by visual 
inspection of the sources rejected in the photometric study of Pa- 
per I. It is worth noting that, besides a dense main-sequence popu- 
lation, this CMD is characterized by a remarkable number of PMS 
stars at faint magnitudes and red colors, as indicated in Paper I. 



investigate the spatial distribution of the identified pre- 
main sequence population, in order to define the limits 
of the main body of the system, and we apply extinc- 
tion measurements. In § |4] we present our calculations 
for new PMS evolutionary models in the observational 
plane. We make use of these models in combination with 
our photometry for the construction of the system IMF 
of LH 95 in §[51 We discuss our findings in §[5] and final 
conclusions are given in § [T] 

2. PHOTOMETRY 

In this study we analyze the photometry presented in 
Paper I, derived from deep observations with the Wide- 
Field Channel (WFC) of ACS on-board HST within the 
GO Program 10566 (PI: D. Gouliermis). Two pointings 
were observed, one centered on the association LH 95 
itself, which we refer to as the "system", and another 
one about 2' to the west on an empty area, typical of 
the local LMC field. We refer to the latter as the "field" . 
For each pointing observations in two photometric bands 
were available: F555W and F814W, roughly equivalent 
to the Johnson V and Cousins I bands respectively. The 
observations and their photometry are described in detail 
in Paper I. We found about 16,000 stars in the system, 
about 2,600 of which fie in the PMS region of the CMD, 
and about 17,000 in the field. 

In the present work we enhance the photometric cata- 
log of stars derived in Paper I by searching for eventually 
missing objects, especially in the low-mass PMS regime. 
Photometry was obtained using the ACS module of the 
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TABLE 1 

Sample from the photometric catalog of all stars found in the 
REGION OF LH 95 with HST/ACS imaging 



R.A. DECL. F555W 0-555 FS14W (7814 

# (J2000.0) {J2000.0) (mag) (mag) (mag) (mag) 



1 


05 37 06.04 


2 


05 37 19.26 


3 


05 36 59.26 


4 


05 37 15.12 


5 


05 37 05.62 


6 


05 36 49.60 


7 


05 37 05.53 


8 


05 37 04.53 


9 


05 37 14.29 


10 


05 36 56.68 


11 


05 36 53.49 


12 


05 36 59.67 


13 


05 37 01.86 


14 


05 36 57.29 


15 


05 36 53.06 



-66 21 37.15 
-66 21 07.88 
-66 21 20.66 
-66 21 44.39 
-66 21 35.39 
-66 23 26.16 
-66 21 59.65 
-66 22 01.02 
-66 22 52.00 
-66 21 10.40 
-66 21 55.55 
-66 22 02.10 
-66 22 35.65 
-66 21 49.03 
-66 22 11.28 



14.839 


0.001 


15.956 


0.001 


15.406 


0.001 


15.503 


0.001 


15.372 


0.001 


16.930 


0.003 


15.632 


0.001 


15.670 


0.001 


15.951 


0.001 


16.828 


0.001 


17.183 


0.002 


16.813 


0.001 


16.752 


0.002 


16.267 


0.001 


16.610 


0.001 



15.066 


0.001 


15.001 


0.001 


15.570 


0.001 


15.494 


0.001 


15.473 


0.001 


14.873 


0.001 


15.740 


0.001 


15.783 


0.002 


15.990 


0.001 


15.102 


0.001 


15.595 


0.001 


14.982 


0.001 


15.663 


0.001 


16.440 


0.001 


16.736 


0.002 



Note. — Magnitudes are given in the Vega system. Units of right ascen- 
sion are hours, minutes, and seconds, and units of decUnation are degrees, 
arcminutes, and arcseconds. The spatial resolution of the ACS/WFC is 
0.05" . 



package DOLPHOT'' (Ver. 1.0). The selection of stars in 
our original photometry was based on the quality param- 
eters estimated for each detected source by the package. 
However, objects which are actual stars may have been 
rejected as spurious detections in crowded regions, espe- 
cially in the presence of a non-uniform background due 
to nebular emission, as well as in the neighborhood of 
very bright sources. Therefore, we considered the catalog 
with all the rejected sources and we performed a visual 
inspection of each of them on the original ACS FITS im- 
ages constructed with Multidrizzle. With this process 
we recovered about 900 additional stars. The updated 
photometric catalog of the association LH 95 includes in 
total 17,245 stars. A sample of this catalog is shown in 
Table [TJ The corresponding CMD is shown in Fig. [TJ 
Typical uncertainties of our photometry as a function of 
the magnitude for both filters are given in Figure (left), 
for both the "system" and the "field" ACS pointings. 

2.1. Completeness 

The completeness of the data is evaluated by arti- 
ficial star experiments with the use of lists of almost 
400,000 artificial stars created with the utility acsfake- 
list of DOLPHOT for each of the observed areas of the 
system and the field. This utility enables to add stars 
of given magnitudes to all the frames simultaneously, 
searching for a detection in none, one, or both photo- 
metric bands. In this way it is possible to derive au- 
tomatically a 2-dimensional completeness function for 
both color {F555W - i^814M^) and magnitude F555W 
or F814W. Considering that at any given faint magni- 
tude in one band the color term covers a range of up to 
2 magnitudes, as shown in the CMD of Figure (TJ allow- 
ing the magnitude in the other band and therefore the 
completeness to vary significantly, this two-dimensional 
approach is necessary for an accurate treatment of the 
completeness of our data. Therefore, we define with our 

* DO LPHOT is an adaptation of the photometry package 
HSTphot l|DolDhirill2000(). The software and its d ocumentation can 
be obtained from http://purcell.as.arizona.edu/dolphot/ 



artificial star technique the 2D completeness as a func- 
tion of (i^555W^ - F8UW) and F8UW, by computing 
the ratio of the detected sources over the total added 
stars in a grid of equally spaced bins 0.25 magnitudes 
wide in {F555W — F814W) and and 0.5 magnitudes in 
i^814W^. 

Such completeness maps have been constructed for 
both the system and the field, enabling, via interpolation, 
an accurate determination of the actual completeness at 
any given point of the CMD. The right panel of Figure [2] 
shows the completeness function in both bands for both 
the system and the field (top and middle respectively). 
The completeness is found to vary with location, and 
therefore different regions in the observed field-of-view 
are characterized by different completeness functions ac- 
cording to their crowding. In the bottom panel of the 
figure we show the completeness as a function of mag- 
nitude for the central region of the association, as it is 
defined in § 13.21 It is remarkably lower than that in the 
entire field-of-view, due to higher crowding of faint stars 
and over-density of the bright ones. 

To avoid further contamination by spurious objects in 
our catalog, which could have been included through ad- 
ditional stars, we trimmed the data set used in this work 
excluding sources with V > 28.5 mag or / > 28 mag. 
At these magnitude limits stars are faint enough to be 
well below any reasonable completeness limit. The com- 
pleteness in the "system" CMD at V=28.5 mag and 1=28 
mag is less than 1%. Consequently, their removal does 
not affect the quality and completeness of our data. 

3. THE PRE-MAIN SEQUENCE POPULATION OF LH 95 

3.1. Field subtraction 

The observed stellar population in the area of the asso- 
ciation LH 95 is naturally contaminated by the general 
field population of LMC. Therefore, in order to study 
the stellar members of the association alone, and in par- 
ticular its PMS stars, the field contamination should be 
removed. In our case, having at our disposal photometry 
in only two bands, the possibility to distinguish if a given 
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Fig. 2. — Left panel: Distribution of the photometric errors, for both the F555W and F814W band for the LH95 system and the LMC 
field. Right panel: Completeness of our photometry i n th e entire observed areas of both the system (top) and the field (middle), as well as 
for the central region of the association LH 95 (see § 13.21 1 . 

the field respectively. Cgys and Cfld are the completeness 
factors, for system and field, measured for exact position 
of the considered star in the CMD, as explained in 12.11 
^ is a normalization factor which takes into account dif- 
ferences in the field stellar density between the two re- 
gions. From repeating measurements in different parts 
of the area of the field and that of the system we found 
that on average the field includes systematically ~ 1.3 
times more main sequence stars than LH 95. Therefore, 
considering that both areas cover equal surfaces, we set 
A = 0.77. According to Eq. [U when P > 1 the star 
is considered to be a field star and it is removed; P = 
means that there are no field stars in the field CMD in the 
elliptical neighborhood of the star and therefore it is con- 
sidered a member of the association. When < P < 1 
the star is randomly kept or rejected as a member star 
with a probability equal to P. 

The result of the first field-subtraction iteration, ap- 
plied to the entire LH 95 field of view is shown in Fig. [31 
where stars flagged to belong to the association are plot- 
ted by black dots. As shown in Paper 1, it is evident 
that the PMS population belongs entirely to LH 95, while 
all lower main sequence stars turn out to belong to the 
LMC field. The upper main sequence (UMS) shows an 
overabundance of stars in the area of the system, demon- 
strating the presence of the corresponding population in 
the association itself. 

3.2. Topography of the PMS Stars 

In this section we study the spatial distribution of the 
PMS stars in the area of LH 95, as an indication of the 



star is a true member or not is limited, but a statistical 
approach can provide accurate results. 

We use a Monte Carlo technique, which considers that 
a star belongs to the cluster or the field with a probabil- 
ity, which depends on the local density in the same part 
of each CMD (of system and field). In order to maxi- 
mize the statistical accuracy of the result, two iterations 
are apphed. First, the field subtraction is applied to the 
entire observed field-of-view centered on the system, en- 
abling to study the spatial distribution of the LH 95 stars 
and to isolate the boundaries of the association. After 
that, the study is limited on this subregion of the frame, 
and another field subtraction is carried out on the central 
region providing an accurate estimate of the membership 
of the sources included. 

Specifically, for the first iteration we consider each star 
in the catalog of the system and we select an elliptical 
region in the CMD of the system centered on the position 
of this star, and the same region in the CMD of the 
field. Each such region is being chosen with semi-axes 
AI = 0.5 mag and A{V — I) = 0.17 mag. This particular 
size has been chosen as to be wide enough to include a 
considerable number of stars, and still small enough to 
achieve a considerably high resolution in the resulting 
CMD. We compute the probability of each star in the 
area of the system to belong to the field as: 



fid 



Cfid 



(1) 



where A'sys and are the numbers of stars included 



in the elliptical region of the CMD of the system and 
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Fig. 3.— F555W - FSliW, F8UW Color-Magnitude Diagram 
(CMD) of the entire observed area around LH 95. The stars which 
are identified as members of the system alone and not of the field 
of the LMC based on the first field-subtraction iteration (see §[3TT]l 
are plotted with thick symbols. The dashed line sets the limits of 
the loci in the CMD of the PMS alone. The spatial distribution of 
these stars has revealed the existence of youn g stellar sub-groups 
within the association (Paper I; see also ^ 13.21 1. 

places where star formation occurs. For this purpose, 
we construct surface density maps of the whole area by 
counting the stars identified as system members after 
applying the Monte Carlo removal of the field contami- 
nation as described in § 13.11 Since we are interested in 
the PMS cluster properties, we considered only the PMS 
population, as selected in the region of the CMD shown 
in Figure [21 

The star counts are performed by dividing the whole 
ACS field-of-view into a grid of 50 x 50 elements, each cor- 
responding to a size of ~ 85 pixels (~ Ipc), and by count- 
ing the (field-subtracted) PMS stars in each of them. 
The selection of the specific grid-element size, which cor- 
responds to the "resolving efficiency" of the detection of 
stellar concentrations, was chosen, after several experi- 
ments, as the most appropriate for revealing the smallest 
concentrations with a physical meaning. 

In order to remove from this map eventual noise, we 
applied a smoothing on the density map with a kernel of 
size ~ 2.5pc ~ 200 pixels. The majority of the PMS stars 
is found to be located within a compact region ~ 1' x 2', 
with a density well above the threshold of 3a, where a 
is the standard deviation of the background surface den- 
sity. The derived two-dimensional density map is shown 
in Figurc|3](left), where the 3a density limit, correspond- 
ing to the first (lower) isopleth, is chosen to define the 
statistically significant concentration of the PMS popu- 
lation of the association LH 95. 

From here on, we refer to the region confined by the 3a 
density isopleth of Figure|4]as the central region of LH 95, 
and we focus our subsequent analysis on the stellar pop- 
ulation included within the boundaries of this region, as 
the most representative of LH 95. In Figure |4] (left) can 
be seen that within the central region there are smaller 
substructures characterized by a higher projected density 
of PMS stars. We select the three most prominent ones, 
isolating circular areas within the central region as shown 
in Figure m (right). We refer to these substructures as 
subcluster A, B and C respectively. The positions and 
sizes of the selected circular regions around the subclus- 



ters, shown in Figure |4] (right), are given in Table [2] with 
those for the whole central region. We measure the to- 
tal mass, Mtot, included in each stellar concentration, 
assuming that all systems follow a mass function sim i- 
lar to the standard Galactic field IMF (" Kroupal 12002^1. 
and extrapolating their stellar content down to 0.08 Mq. 
We then estimate the stellar density, g, and the disrup- 
tion time, td, of each c luster following the method by 
iGouliermis et al. (120021 their section 6). The latter is 



given as (lSpitzedfl95ak 



td = 1.9 X lO^g 

\pc3 



years 



(2) 



Moreover, the dynamical status of a stellar system is de- 
fined by two additional ti me-scales, the crossing and the 
two-body relaxation time ()Kroupall2008[ ). which are given 
as: ^ ^ 

t„ = and trolax = 0.1. — tcr (3) 

a in A* 

respectively. The three-dimensional velocity dispersion 
of the stars in the cluster, cr, is given as 



GMto 



(4) 



where e is the star formation efficiency (SFE) and rh 
the half-mass radius of the cluster. In order to make a 
rough estimation of the aforementioned time-scales for 
the subclusters and the whole central region of LH 95, 
we apply the formulas of Eq. [3l 

The SFE in several nearby Galactic gas-embedded 
clust ers has been found t o range typicall y from 10% t o 
30% (|Lada fc Lada ll2003l) or 20% to 40% (|Kroupall2008fl . 
value which increases with time while the gas is removed. 
However, the average age of the system is 4 Myr (see !JS]), 
slightly greater than t he typical time nec essary to re- 
move most of the gas (jLada fc Lada 1 120031 ) . Therefore, 
and taking into account the very low optical extinction 
(see i )3.2p . LH 95 should be considered as an emerging 
cluster, meaning that it should be at the process of sep- 
ar ation from the pa r ental cloud. Furthermore, according 
to IWilking fc Ladal (|1983f ). a high SFE is required for a 
bound young cluster to emerge from its parental cloud. 
As a consequence, we consider a value of e ~ 0.4 as more 
adequate for our system. We derive, thus, values of a 
of the order of 2 km s~^ for our objects and we pro- 
vide the additional estimated structural parameters for 
each system also in Table [2j We find that the crossing 
time for the whole region is greater than the age of the 
system, and this is quite consistent with the observed 
sub-clustering of the system. 

From a comparison between the characteristics of the 
subclusters and the central region, as they are given in 
Table [2l one can see that there is a significant fraction 
of mass outside the immediate regions of the subclus- 
ters. This can also be seen in the maps of Figure 21 
where a prominent population of PMS stars is easily dis- 
tinguished away from the subclusters, spread in the area 
among them within the central region. We estimate that 
the fraction of distributed stellar mass over the total cor- 
responds to ^ 40%, with the remaining 60% being clus- 
tered. This fraction of distributed PMS stars in LH 95 
is consistent with statistical analyses of several Galactic 
nearby star-forming regions that predict ^ 60% (see 
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Fig. 4. — Left panel: Isodensity contour map of the area of the system constructed from star counts of the PMS stars alone. Isopleths 
are drawn in steps of la, where a corresponds to the standard deviation of the background density, starting from the 3cr level. This level 
defines the limits chosen for the selection of the stars used for the construction of the CMD of the central part of the area of the system 
shown in Figure[5]{left and right panels). This map demonstrates that, as we found in Paper 1, PMS stars are concentrated in few compact 
subgroups within the association. These density peaks in the PMS stars coincide with those of the Upper Main Sequence (UMS) stars 
(Paper I), suggesting that these subgroups contain fully populated IMFs. Units are in pixels on the ACS drizzled frame. Right panel: 
Positions of the three most prominent concentrations (subclusters) of PMS stars, which appear as the densest substructures of the main 
part of the association LH 95. 



TABLE 2 

Characteristics of the Central Region (LH 95) and its three subclusters 





RA 


DEC 


Size 


rh 


Aftot 






^cr 


^rclax 




(J2000.0) 


(J2000.0) 


(pc) 


(pc) 


(103 Mq) 


(Mq pc-3) 


(100 Myr) 


(Myr) 


(Myr) 


LH 95 


5 37 06.50 


-66 22 03.50 


33.5 


6.5 


2.37 


0.47 


0.9 


6.4 


433.0 


A 


5 37 15.25 


-66 21 40.87 


8.7 


1.7 


0.33 


1.28 


2.4 


2.2 


26.9 


B 


5 36 59.78 


-66 21 36.88 


4.7 


0.8 


0.19 


2.70 


5.1 


0.9 


6.7 


C 


5 37 03.96 


-66 22 09.73 


12.9 


3.0 


0.98 


1.11 


2.1 


3.1 


95.5 



review bv lAllen et al.l [2"007[ ) . However, it is not clear if 
this distributed population is the result of the star for- 
mation or a merging process. Detailed simulations do 
predi ct that subclusters m ay merge to form a larger one 
(e.g. iFellhauer et"al] l2006f) . but our data do not allow 
us to verify if this takes place in LH 95. In any case, 
the stellar density and disruption time of all considered 
systems, including the whole central region, show that 
all are rather compact and none of them is under disrup- 
tion. Subclusters A and B, located at the east and west of 
the central region, are particularly compact and roughly 
spherical shaped, while subcluster C, close to the south- 
ern limit of the central region, presents a more elongated 
and extended distribution of PMS stars. This subcluster 
corresponds to the visibly prominent main part of the 
association (see Figure 1 in Paper I). 

3.3. LH 95 Central region 

For the detailed study of the PMS population in LH 95, 
while minimizing the contamination by other stars, we 
focus on the central region of the system. Specifically, 
we consider the entire population included in the central 
region and we apply, for a second time, the Monte Carlo 



technique for field subtraction, taking into account the 
ratio between the surface of the observed LMC field and 
that of the central region. 

It should be noted that this second application of the 
field subtraction technique to the central region is not 
equivalent of simply isolating the stars, which are found 
to be located within the boundaries of the region after 
the first iteration of the Monte Carlo subtraction. 

Naturally, the central region has a higher density 
of system stars (specifically PMS stars), and therefore 
a lower relative contamination by field stars. Conse- 
quently, the actual fraction of stars which are marked to 
be field stars is lower, reducing the probability of (statis- 
tically probable) wrong membership assignments. The 
number of field stars included in the elliptical region of 
the CMD TVfid is computed based on the photometry in 
the entire field and scaled according to the area occu- 
pied by the LH 95 central region, which we found to be 
~ 12.6% of the total ACS field-of-view. 

The result is shown in Figure [5l where we plot the 
CMDs of the whole population included within the 
boundaries of the central region (left panel), the one of 
the whole area of the field, sampled randomly and down- 
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Fig. 5. — Left: CMD of all stars observed in the central part of the area of the system. This part was selected as the most representative 
of the association based on star counts (§ 3.2; Figure|4]l. Middle: CMD of the corresponding part of the area of the field, which was used 
for the decontamination of the CMD of the system from the contribution of the field stellar population. Right: CMD of central area of the 
system after the field subtraction was applied with the use of an advanced Monte Carlo method (§ 3.1). This CMD highlights the UMS 
and PMS stars as the prominent stellar populations of the association LH 95, demonstrating its youthfulness. 



scaled according to the area coverage of the central re- 
gion (central panel), and the corresponding CMD of the 
central region after the field contribution has been re- 
moved (right panel). From this figure it is evident that 
the observed stellar population of the low main sequence 
(LMS) can be entirely considered as field population. On 
the other hand the association has a prominent upper 
main sequence (UMS) component, and all pre-main se- 
quence stars, easily distinguished by redder colors, are 
solely present in LH 95 alone. Furthermore, the second 
application of our field subtraction technique, limited to 
the central region, allows us to reduce considerably the 
contamination by objects not located in the PMS and 
UMS parts of the CMD, which are still present as resid- 
uals from the first application of the method on the entire 
field-of-view (Figure [3]). 

In the CMDs of Figure \5\ can be noted that the field 
population has a tighter sequence than the PMS popula- 
tion, which may indicate an age-spread among the PMS 
stars of LH 95 (see also Figure ^ . Indeed, the loca- 
tions of low-mass PMS stars in the CMD of young stel- 
lar systems of our Galaxy often show a widening, which 
may be evidence for an age-spread in the system (e.g. 
iPalla fc Stahleill2000f ). However, several characteristics 
of these PMS stars (being T Tauri stars), such as variabil- 
ity and circumstellar extinction, can cause considerable 
deviations of their positions in the CMD (jSherrv et al.l 
[2004), which may be misinterpreted as an age-spread. 
As a consequence, detailed simulations of the character- 



istics of PMS stars are required to quantify the effect of 
these characteristics on their positio ns in the CMD and 
to co nclude on any true age-spread (jHennekemper et al.l 
I2008f) with the use of population-synthesis techniques. 

In the following sections we focus our study on the 
isolated, field subtracted population in the LH 95 central 
region, shown in Figure [5] (right). 

3.4. Interstellar Extinction 

With photometry available only in two bands it is not 
possible in general to measure the visual extinction Ay 
independently for every star of the system. It is, however, 
possible to perform a statistical approximation obtaining 
the average value of reddening as well as its distribution. 
For this purpose we consider the UMS stars (/ < 20.0 
mag) in the central part of LH 95, after performing field 
subtraction, and we compare the position observed in 
the CMD to that expected according to the Padova grid 
of evolutionary mod els in the ACS photometric system 
(|Girardi et al.H200l . Considering that LH 95 is a very 
young association (r ^ 10 Myr) we assume that all UMS 
stars have a age equal to the youngest available isochrone 
in the Padova models, logr — 6.6 for a metallicity of 
Z = 0.008 (typical value for the LMC; iKontizas etall 
1199 ^ and distance modulus m — M = 18.41 mag, com- 
patible with several distance m easurements for the L MC, 
available in the literature (e.g. lAlcock et al.l (|2004f )). 

For every star we compute the intersection between the 
reddening vector and the isochrone, and we obtain, thus. 
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Fig. 6. — Determination of the reddening in the region of LH 95. Left: The field subtracted stellar population within the central area 
of LH 95 with the isochrone model for logjQ r = 6.6 and the corresponding individual reddening vectors of the upper main sequence stars 
overlayed. Square symbols represent objects with {V — I) > 0.3 mag, which may be considered as candidate Herbig Ae/Be stars, due to 
their highly reddened colors. Middle: Distributions of the measured reddening for all UMS stars (filled histogram) and excluding candidate 
Herbig Ae/Be (darker histogram). Right: UMS CMD showing the Girardi et al. (2002) isochrone for E{F555W - F8UW) = and with 
the mean reddening from the distributions shown in the central panel applied. 



the reddening value in terms of E{F555W — FSliW). 
Since the ACS photometric system is shghtly different 
from the standard Johnson-Cousins system, we compute 
the exact reddening parameters for our ACS ba nds us- 
ing the typical galactic extinction law of Cardelli et al.l 
(|1989[ ). parameterized by a value of Ry = Av/E{B — 
V) = 3.1. We consider a s et of template spectra taken 
from the NextGen catalog (|Hauschildt et al.l[l999f ). and 
we apply the extinction curve with an arbitrary Ay 
which we choose to be equal to 1, and we measure 
the extinction Ap^^^w and Ap^nw by means of inte- 
gration of the original and reddened spectra within the 
ACS filter profiles. In this way we compute Rf555W = 
Af555w/eIf555W - F8UW) ~ 2.18 and Rfsuw = 
Af8uw/E{F555W - F8UW) ~ 1.18. 

The derived reddening distribution, shown in Figure [H] 
(middle), appears to have a tail extended to higher val- 
ues (filled histogram). This tail appears due to a num- 
ber of UMS objects, which show a high color excesses. 
Taking into account all objects, including the ones with 
high color excess, the average reddening has a value 
E{F555W - i^814M^) ~ 0.34 mag, corresponding to an 
optical extinction of Ay = 1.56. However, these red ob- 
jects are being suggested to be Herbig Ae/Be (HAeBe) 
stars (|Gouliermis et al1l2002l ). Such stars of intermedi- 
ate mass, aged between 1 and 10 Myr, occupy the faint 
end of the UMS with 17.5 < F8UW < 20 mag. Being 
in their pre-main sequence phase, they have an intrinsic 
color higher than that of MS stars in the same luminos- 
ity range. Consequently, if indeed these are not MS but 
HAeBe stars, which is possible due to the youthfulness of 
the system, then the derived reddening would be biased 
towards higher values. 

We estimate the mean reddening without considering 
these stars (darker histogram in Figure [6] - middle) and 
we derived a lower mean value of E {F555W — FSIAW) ~ 
0.275 mag, corresponding to an optical extinction of 
Ay = 0.6 mag. In Figure [6] (right) both values are ap- 
plied to the 4 Myr isochrone model (the youngest avail- 
able in the grid of models by the Padova Group). From 
this figure it is evident that the difference between the 
two values is not dramatic, producing more or less the 
same shift of the model in the CMD. As a consequence. 



and in order to avoid any biases due to the presence of 
young stellar objects in LH 95, we consider an optical 
reddening of E{F555W - i^814W^) ~ 0.3 mag as the 
most representative for the region. We repeated the con- 
struction of the reddening distributions shown in Figure[6] 
(middle) and the determination of the mean reddening 
values using different isochrones, spanning an age range 
between 4 and 7 Myr and assuming different distances 
varying by 0.2 mag around the assumed distance mod- 
ulus. We found that the mean reddening remains prac- 
tically unchanged within 0.02 magnitudes, a difference 
small enough to confirm the accuracy of our estimation. 

4. NEW SET OF OBSERVATIONAL PMS MODELS 

In the last decades several theoretical models for the 
pre-main sequence stellar evolution have been computed. 
Some of the most popular, for PMS masses above the 
hydro ge n burning limit , are t h e models of ISiess et al 



(f2000), 'P alla k Sta bler (ToM), 'D'Anto na fc Mazzitelli 
(1994, 1993) and [Swcnson et al. (1994). These models 
are generally expressed in terms of physical quantities, 
such as the effective temperature and the total bolomet- 
ric luminosity, describing the evolution of stars in the 
HR diagram. However, one of the primary aims of stellar 
evolution theory is the explanation of the observed pho- 
tometric data of stars in order to extract their masses and 
ages from their magnitudes and colors, and therefore a 
conversion between physical and observable quantities of 
the models i s req uired. Such a conversion is specifically 
described by Sfess et al. (2000), who include the transfor- 
mation between Tos and L to colors and magnitudes in 
the U BVRI Cousins system and JHKL infrared bands 
for their models. These authors use simple relations T^s 
versus color, as well as bolometric corrections from either 
Siess et al. (1997) or Kcnyon & Hartmann (1995). These 
conversion tables, derived from observations of stellar 
clusters, are valid for solar metallicity dwarf stars, but 
the discrepancies can be significant when one deals with 
populations with different stellar parameters. In partic- 
ular, the age of a star is related to the stellar radius and 
therefore to the surface gravity, which introduces differ- 
ences in the spectral behavior, and therefore a popula- 
tion of different ages or metallicities may require differ- 



The IMF of LH 95 in the LMC down to the Sub-Solar Regime 



9 



lo'r 



q(q) 



4 

[cqsl 



Fig. 7. — Distribution in the T^g, log (g) space of tlie syn- 
thetic spectra collected for our conversion of theoretical evolution- 
ary models into the observable plane. The same spectra are used 
for metallicities of [M/H] = 0, —0.5, —1, —1.5. Squares are spec- 
tra taken from the Next Gen models of iHau schildt ctal. (1999), 
while crosses are those of lKuruc3 II1979I . [l993f) (see ^ 14.11 1. 

ent conversion relations. This issue can be important for 
e.g. cold M-type stars, for which the broad molecular 
absorption bands dominate the optical spectra, affecting 
integrated colors and color corrections. 

A more thorough method to analyze this issue is to 
make use of synthetic atmosphere models, performing 
hoto metry directly on synthetic spectra. Girardi et al.l 
20021 ) applied this method for their evolutionary models 
for evolved populations, which were converted in sev- 
eral photometric systems using a grid of atmosphere 
models described by three parameters {[M/H], log (Tes), 
log (5)). This method was also ap plied for stars o f sub- 
solar mass and brown dwarfs by iBaraffe et al.l (|1998l 
|2001|) . For low- and intermediate-mass stars, however, 
there are no similar conversions of evolutionary models 
into the observable plane for various photometric systems 
available. As a consequence, and in order to have an ac- 
curate transformation of the PMS evolutionary models of 
Siess et al. (2000) into the observational plane, we apply 
a conversion o f these models follow ing an approach sim- 
ilar to that of iGirardi et al.l (|2002f ). We construct, thus, 
observational models (both tracks and isochrones) for 
four assumed metallicities (Z = 0.02, 0.01, 0.008, 0.004) 
and several photometric systems (specifically Johnson, 
Cousins, HST WFPC2, HST ACS, HST STIS^ HST NIC- 
MOS, Galex, Stromgren, SDSS, Tycho and 2MASS). 

4.1. The Method 

We apply our conversion method first by collecting a 
grid of atmospheric models in the T^s, log{g), [M/H] 
space. These models are shown in Fig. [71 We selected 
spectra from NextGen (|Hauschildt et aLri999|) for low- 
and medium-temperature stars. The selected grid covers 



temperatures 2600 K < Toff < 8200 K in steps of 200 K 
and surface gravities 2.0 < log {g) < 5.5 in cgs units, 
spaced in 0.5 dex. For higher ternperatu res we used the 
grid of models by iKurucd ()1979l . I1993D , with gravities 
log (g) between 3.0 and 5.0. This grid covers entirely the 
range of stellar parameters in the evolutionary models of 
Siess et al. (2000). The considered metallicities cover 
the values [M/H] = 0.0, -0.5, -1.0,-1.5, representa- 
tive of stellar populations in the Magellanic Clouds, as 
well as in many regions of the Galaxy. We re-sampled 
the NextGen models, the resolution of which is much 
higher than required for our purposes of broad-band syn- 
thetic photometry, to the wavelength values of the Ku- 
rucz models, and we normalized them to the same units 
(erg cm~^ s~^ A~^). 

The absolute magnitude in a photometric band S\ of 
a star with a spectral energy distribution Fx , and stellar 
radius R is given by 



Ms. 



-2.5 log 



R 
lOpc 



AFaS'aIO 



dA, 



-ZP; 



(5) 

where is a reference spectrum that gives a known 
apparent magnitude ZPs.; in the Vegamag photomet- 
ric system, which uses the flux of a Lyr as reference, 
fx = Fx^vEGA and the zero-points ZP are close to zero^. 

We assume an extinction Ax = 0, since we are inter- 
ested in constructing a general observable grid of evolu- 
tionary models for unreddened objects. Eq. ([5]) can be 
written then as: 



Ms 



where 



B{Fx,Sx 



43.2337- 5 log 



R^ 

i?0 



+ BiFx,Sx) (6) 



BiFx,Sx) ^~2.5log 



XFxSxdX 



XF: 



X.VEGAJX 



SxdX 



ZPs. 



ij) 

The latter term can be directly calculated having the 
synthetic spectrum of any considered star, a calibrated 
Vega spectrum, the band profile of the photometric filter 
considered and the eventual zero-point ZP. For the com- 
putation of B we used the Gensynphot code i ncluded 
in the Chorizos package (|Maiz-Apellanid[200l ) . Given 
a set of stellar parameters (Teff, log (5), and [M/H]) for 
a sample of points (stars) on an evolutionary track or an 
isochrone, the corresponding spectra are computed by 
interpolation on the Cartesian grid of synthetic models 
of stellar atmospheres. Then, the integral at the numer- 
ator of Eq. ([7]) is computed. The final normalization 
is performed using a recent reference spectrum of Vega 
(|Bohlinll2007D with the co rrected zero-points as reported 
bv lMaiz-ApellanizI ()2007[) . 

^ Although the historical Vegamag photometric system was sup- 
posed to have zero-points defined so that to impose the apparent 
magnitude of Vega equal to zero, improvements in the modeling 
of the spectrum of Vega, in its measure and in the definition of 
standard throughputs, have introduced small corrections. 
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Fig. 8. — Comparison of the evolutionary models for PMS stars de veloped in this work with the ones bv lSiess et"an ll2000fl . Left: 
Color-Color Diagram of isochrones from the models of 'Sies s et all ll2000h for the ZAMS and ages of 0.1, 0.3, 1, 3, 10, 30 Myr for Z = 0.02. 
The dotted line refers to the original Siess et al. (2000) isochrones, in which a unique temperature — color relation is used, while the solid 
lines are the results of our analysis. The differences between the two treatments are mostly evident at larger colors, at about T ~ 4000ii', 
where M-type optical spectra are dominated by oxide bands, and a re strongly dependent on age due to variations in the surface gravity 
of the stars. Right: Comparison between the original isochrones of |Siess_et_aI.' (2000) for Z = 0.02 and our treatment in the {V — I), I 
diagrams. For medium masses (top left region) the differences are negligible, while for l ow-mas s , col d stars the shifts are important. In 
particular it is evident that our transformations extend to larger colors than those of Siess et al.l 1I2OOO1 '). 



Following this method we converted al l the orig- 
inal evolutionary models of ISiess et all (|2000f ) into 
absolute magnitudes, for the metallicities oi Z ^ 
0.02, 0.01, 0.008 and 0.004 corresponding to [M/H] ~ 
0.0, —0.3, —0.4, —0.7 respectively. Considering that 
the original Siess et al. models were computed only 
for Z = 0.01 and 0.02, we extrapolated the theoretical 
isochrones computed for Z = 0.01 to lower metallicities 
with the use of the atmosphere models computed for such 
metallicities. Naturally, a more detailed treatment of 
low-metallicity tracks and isochrones requires both evo- 
lutionary and atmospheric models to be available at the 
same considered metallicity. However, for the present 
analysis in the LMC, for which Z ^ 0.008, the uncer- 
tainties are quite low given the small difference between 
the metallicity of the theoretical isochrones (0.01) and 
that of the atmosphere models (0.008). 

4.2. Comparison with previous conversions 

With our construction of the observational plane for 
PMS evolutionary models, we verified that a single 
temperature-color relation is not sufRcicnt for stars with 
different physical parameters, especially surface grav- 
ity. In Figure [5] (left) we present a comparison between 
the {B — V), {V — 1) color-color diagrams for differ- 
ent isochrones, as they were derived from the method 
of ISiess et all (2000) (red dotted lines) and from our 
method with the use of atmosphere models (blue contin- 
u ous lines). In this f igure it is shown that the isochrones 
of ISiess et all ([20001 ) are distributed along one single se- 
quence, due to the use of monolithic temperature-color 
relations, which are independent of the surface gravity. 
However, our observational isochrones show a spread for 
high color terms (low Tcff), because differences in ages 
imply differences in the surface gravities, which produce 



variations in the integrated colors according to the at- 
mosphere models we used. The major differences arise 
for M-type stars, in which the optical spectra are domi- 
nated by strong oxide absorption bands. Figure [8] (right 
panel) shows these differences on the CMD. As shown 
in this fi gure, the differe nces between our method and 
that of ISiess et al.l (|2000( ) are mostly evident at the red 
end of the models, where our computations predict red- 
der colors, and for the youngest models, where the stellar 
surface gravities are lower than those for main sequence 
stars. In Figure [8] the example oi Z — 0.02 is shown. 

5. THE INITIAL MASS FUNCTION 

The mass function (MF) is defined as the number 
distribution of stars as a function of mass. While a 
general observational approach measures the so called 
Present Day Mass Function (PDMF), a cornerstone for 
understanding how stars form is the Initial Mass Func- 
tion (IMF), which is the mass distribution according to 
which stars are born. For a pre-main sequence popula- 
tion younger than the time required by the most massive 
stars to disperse it, or to evolve into post-MS evolution- 
ary phases, the observed mass function coincides with 
the IMF. 

Generally, the IMF is parameterized as follows: 



^(M) dM cx M~(i^ 



(8) 



name ly, approximating it with a power-law (e.g. lSalpeteil 
11955( 1 or with a series of power-l aws, with exponents 
chang ing in different mass ranges l|ScaIol[l98fit iKroupal 
120021 ). In this section we describe the "average" IMF de- 
rived for the stellar populations of the central region of 
LH 95, as it has been defined in S 13.21 

The F555W-F8UW, F8UW CMD of Figure [5] (right 
panel) is shown again in Figure [9l with our observable 
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PMS evolutionary tracks (left panel) and three indica- 
tive PMS isochrones (right panel), computed in §|4l over- 
layed. In the right panel of the figure the ZAMS from our 
grid of models is also plotted along with the youngest MS 
isochrone from the Padova grid of evolutionary models 
(jGirardi et al.ll200^ corresponding to logr ~ 6.6. Both 
isochrones fit to each other very well demonstrating that 
the bright MS stars of LH 95 have ages of ^ 4 Myr. As 
far as the PMS stars of LH 95 concerns, their positions in 
the CMD corresponds to an age-spread that covers ages 
between ~ 1.5 and 10 Myr according to our PMS models. 
A simple statistics on the distribution of the ages of these 
PMS stars, based on their CMD positions, results to a 
normal distribution of ages which peaks at about 4 Myr. 
However, if this age-spread is true or not can be the sub- 
ject of debate, as we discuss in § 13.31 Therefore, based 
on the PMS population of LH 95 we define an indicative 
age for LH 95 of logr = 6.6 ± 0.4. Corresponding PMS 
isochrones for ages of 1.5, 4 and 10 Myr from our grid of 
observable PMS models are overlayed on the CMD in the 
right panel of Figure [5] A comparison of this age, which 
also represents the individual subclusters of LH 95, with 
the corresponding dynamical crossing time, estimated in 
§ 13.21 shows that all sub-clusters, as well as the whole 
central region have an age comparable to their crossing 
times, indicating that these systems are under 'mixing' 
process if not already 'mixed' and therefore we probably 
do not observe them close to their initial morphology. 

5.1. Construction of the IMF 

Based on our new set of PMS models computed above 
for the ACS photometric bands and the average LMC 
metallicity {Z ~ 0.008), we assigned a value of mass to 
all the stars identified in the region of LH 95, by com- 
paring our observed data with the evolutionary tracks 
of these models. For this comparison we assumed a dis- 
tance modulus of m — Af = 18.41 mag, and an average 
reddening of E{F555W - i^814W^) = 0.275 mag, as it 
is derived in § 13.41 According to Figure [51 our photom- 
etry reveals PMS stars with masses as low as 0.2 M0. 
Our sample, though, as we discuss later, is complete for 
masses down to ~ 0.43 Mq. For the PMS stars we assign 
a mass measurement to each observed star according to 
their positions in the CMD by interpolating between the 
PMS evolutionary tracks shown in Figure [9l 

The position of the brightest main sequence stars is 
not covered by our PMS evolutionary models, limited 
for masses below 7Mq. Therefore, we derive the masses 
of these stars directly from the Padova 4 Myr isochrone 
- already used for the determination of the average red- 
dening in § 13.41 - shifted in the direction of the redden- 
ing vector to reach individually each one of these stars. 
The good agreement between this Padova isochrone and 
the corrected ZAMS from the Siess models is evident 
from Figure [9l While we de-redden the most massive 
MS stars according to their individual position in the 
CMD, their mass estimation is based solely on pho- 
tometry. A more accurate estimation of the mass for 
massive stars certainly requires spectroscopy, because 
for such hot stars t he broad-band colors are degenerate 
(jMassev et al.l[l"995l ). and naturally this effect may influ- 
ence the mass estimation for the brightest MS stars in 
our sample and the corresponding IMF. We discuss the 
influence of the use of photometry alone to the slope of 



the derived IMF in fJSH 

With the masses of both the MS and the PMS stars in 
our sample available the construction of their IMF is in 
principle straightforward, by binning the stars according 
to their masses and fitting the derived distribution. How- 
ever, there are considerations, which affect the construc- 
tion and interpretation of the IMF, and which should be 
taken into account. We discuss these issues below. 

(i) An important issue to be considered is the decon- 
tamination process of the stellar sample in the re- 
gion of the association from the field populations. 
The Monte Carlo method for the field subtraction, 
as applied in § 13. 3i is based on probabilities and 
therefore it is a stochastic process. As a conse- 
quence, the number of stellar members of the sys- 
tem, which are counted in an arbitrary mass bin is 
not determined in an unambiguous manner. More- 
over, the statistical significance of this number is 
not determined by the fraction of stars considered 
as system members over the total number of stars, 
but it depends also on the probability that all the 
observed stars in the bin are system members. This 
dependence affects specifically the MS populations 
with 19 i$ mF8i4w/mag ^ 22, which as it is 
shown in Figure [5] are well mixed with field popu- 
lations, implying a higher uncertainty in the deter- 
mination of the fraction of system members. In the 
faint part of the CMD, system and field populations 
are in general well distinguished, and therefore the 
uncertainties in the numbers of PMS stars that 
are true system-members are insignificant. How- 
ever, there is an important possibility that a frac- 
tion of the binary sequence, evident in the CMD 
for the LMC field population (see Figures [5] and 
HI, contaminates the selected sample of PMS star- 
members of the system. We consider this bias and 
report on its effect in § 15.21 

(ii) The fitting process of the IMF. Typical linear re- 
gression methods, which are used for obtaining a 
functional form of an observed IMF, such as the 
X^-method, consider the presence of measurement 
errors and are based on the assumption that the un- 
certainty associated to each point follows a Gaus- 
sian distribution. However, in this case the uncer- 
tainty in the number of counts within a mass bin is 
the overall effect of both the Poissonian error that 
naturally comes from the counting process, and the 
uncertainty that arises from the field subtraction, 
as discussed above. Although the first can be well 
approximated with a Gaussian for large n, the lat- 
ter depends on the positions of stars and tracks 
in the CMDs. Moreover, in the high mass regime, 
where the number statistics is poor, the Poissonian 
error cannot be treated as Gaussian, requiring the 
asymmetry of the distribution to be taken into ac- 
count. 

(iii) The completeness-correction process. In the low- 
mass regime corrections for the incompleteness of 
the stellar sample are required to estimate the ac- 
tual number of stars. However, given the width 
of the low-mass evolutionary tracks in the CMD 
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Fig. 9. — The CMD of the cluster area of LH 95 with PMS evolutionary tracks (left) an d iso chrones (right) overlayed for a distance 
modulus m — M = 18.41 mag and a mean reddening of E(Fb5bW — FfilAW) = 0.275 mag (see i) l3.4l l. The PMS tracks and the corresponding 
PMS isochrones and ZAMS are computed by us for the ACS filter system, and for the metallicity of the LMC (§|4jl. The 4 Myr isochrone 
for evolved MS stars is from the Padova grid of evolutionary models I Girardi et al.|.2d02, ). The comparison of this isochrone with our ZAMS 
demonstrates that it practically corresponds to the youngest MS populations. The selected PMS isochrones overlayed in the right panel 
have ages of 1.5, 4 and 10 Myr. The central, 4Myr, corresponds to the the mean age of the observed PMS stars, where the other two 
delimit the Icr width of the broadening of these stars in the CMD. 



shown in Figure [H and the 2-dimensional varia- 
tion of the completeness itself in the CMD (in both 
magnitudes and colors; see ij |2.ip . this correction is 
not unique for all the stars counted in a considered 
mass bin. 



We construct the system IMF of the association LH 95, 
with a special care in addressing the above issues, as 
follows. We first determine a grid of logarithmic mass 
bins, for which we choose variable width in log(Af), 
that increases from 0.04 dex at low masses to 0.25 dex 
for the higher masses. The counting of the stars is 
made in variable-sized bins because such bins yield very 
small biases, which are only weakly dependent on the 
number of stars, in contrast to uniformly binned data 
|Maiz-Apellaniz fc Ubedall2005D . For each z*'' bin a num- 
ber of stars Ni is derived. In order to achieve a statisti- 
cally correct sampling of points for the constructed IMF, 
we simulate for each bin the probability distribution ex- 
pected for every Ni with a sample of 1000 points ran- 
domly selected from a Poisson distribution with mean 
N,. 

We then consider the position in the CMD of all the 
stars in each bin, and we compute for each star its com- 
pleteness C (the ratio between the number of sources de- 
tected and the total number of sources) as a function of 
both magnitude and color according to the completeness 
measurements described in § 12.11 for the central region 
of LH 95. Let us consider the j'^ star counted in the i*^ 
bin. The star is located in a specific position of the CMD, 
where completeness is Cj, meaning that the total num- 
ber number of stars corresponding to j is simply l/Cj. 
Consequently, for the total number of stars counted in 



the bin, Ni, the completeness correction is: 



N, 



(9) 



which is a quantity > 1. 
We multiply our Poissonian sample associated to each 
bin by Fi, and we divide by the width of the bin in 
order to express numbers in units of number of stars per 
solar mass^. In order to take into account the additional 
error introduced by the stochasticity of the Monte Carlo 
field-subtraction 13. 3[) . we iterate the technique of field 
removal 200 times, each time repeating the procedure 
explained in the previous paragraphs, namely, counting 
the Ni number of stars in each i*'' bin, sampling a Poisson 
random population associated to it, and computing and 
applying the completeness correction. 

The outcome of this approach is that for each mass- 
bin, instead of a single point-error combination, we have 
a population of 200,000 random IMF points, derived ac- 
cording to the true statistical distribution of the expected 
number of stars in the bin. Clearly, the mean and the 
standard deviation of each of these distributions repre- 
sent the IMF data points and error bars. The main ad- 
vantage of this technique is that we override the non- 
Gaussian behavior of the mass-bin error bars. Specifi- 
cally, now for every mass-bin we do not have one point 

^ It is worth-noting that the completeness correction should be 
applied after simulating the Poisson distribution, and not before. 
For example, if one star of 25% completeness is counted into the 
i"^ bin, the true number of sources associated to the Ni sample is 
4, with an uncertainty of \/T * 4 = 4 and not \/i = 2. 
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Fig. 10. — The system IMF plot of LH 95, which demonstrates 
our method for treating the distribution of the expected number 
of stars per mass-bin according to uncertainties due to counting, 
incompleteness and field subtraction. Each mass-bin is represented 
by the distribution along the j/— axis of a set of 200,000 randomly 
selected points, which map the true statistical uncertainties of star 
counts. An additional Gaussian spread with cr = 1/8 of the bin 
width along the a;— axis has been artificially added for each bin to 
better highlight the behavior of the individual distributions. The 
asymmetry of the statistical distributions especially for the high- 
mass end of the IMF is demonstrated by the drop-like shape of the 
distributions. The very-low mass end (light-gray points) appears to 
drop due to the limitation of our completeness correction technique 
in CMD regions where no stars are detected. Units of the IMF are 
number of stars per solar mass per pc'^ . 

associated with an error that cannot be treated by the 
regression theory, but a numerous set of points with no 
uncertainties associated to them, and consequently with 
arbitrary equal weights for all mass-bins. Naturally, this 
allows a more accurate fitting of the mass function. 

In Figure [TUl we visualize the concept of this approach. 
For every mass-bin the set of random points, represent- 
ing the distribution of the expected S,{M) is drawn with 
tiny colored dots. In this plot an artificial spread along 
the X-axis (a random Gaussian shift with a = 1/8 of 
the bin width) is introduced to better highlight the dis- 
tributions along the y— axis. The black dots with error 
bars show the average and standard deviation of each 
sample of points. The derived numbers are reliable for 
M <; 0.35 Mq, whereas for the lower masses (first five 
mass bins from the left) the completeness corrections are 
not sufficient. This limitation is evident from the CMD 
of Figure [H where it can be seen that for the evolution- 
ary tracks for very low masses {M ^ 0.35 Mq), the 
completeness is extremely low, so that no stars can be 
detected, limiting the ability to recover the true number 
distributions for such masses. 

5.2. IMF model fitting 

A commonly used procedure to derive the functional 
form of an observed IMF is to fit a predefined function 
to the individual numbers of stars originally counted per 
mass-bin taking into account the corresponding count- 
ing errors. However, in our case, we have simulated the 
statistical distribution of ^(M) per mass-bin taking into 



13 

account the additional errors due to incompleteness and 
field subtraction. Within our statistical approach we per- 
form 200 repetitions of the Monte Carlo field subtraction 
method and for each of them we map the Poissonian 
statistics 1000 times. The outcome is a non-analytical 
probability distribution for each mass-bin that represents 
the true expected statistical uncertainties due to the is- 
sues discussed in This process provided a large 
sample of 200,000 points per mass-bin, and naturally the 
best-fitting function can be derived more accurately with 
the use of these data, rather than the original 14 bins. 
The only assumption that we make is that the functional 
form of the IMF of LH 95 can be correctly represented 
by power-laws, according to the parameterization given 
in Eq. ([8]), and therefore we search for the best multi- 
ple power-law which may reproduce our data. As seen 
in Figure [10] the IMF seems to be more shallow for the 
low masses (M < 1 Mq) than for intermediate and high 
masses, implying that it could be reproduced by a two- 
phase power-law. We investigate, however, whether a 
two-phase power-law is sufficie nt or a three- phase power- 
law, such as that proposed by IScalol (|1998D . is required 
to represent our IMF. 

In order to find the best fit of two- and three- 
phase power-laws to our IMF we applied a Levenberg- 
Marqu ard non-linear l e ast square minimization tech- 
nique (|Levenberd 119441 : iMarguardtl [19631 . For this fit 
we consider only the mass-bins of the IMF for which the 
completeness correction allowed us to reconstruct the ac- 
tual number of stars, meaning for log(M/M0) > —0.5. 
The slopes of the power-laws, as well as the position of 
the break points along the abscissa are the free param- 
eters of our fit algorithm. All 200,000 points of each 
mass-bin are considered of equal weight. However, since 
a two-phase power-law is a sub-case of the three-phase 
power-law, the sum-of-squares of the latter is always 
lower. Therefore, the results of the fit algorithm for each 
case are normalized in a rigorous manner by performing 
the so-called statistical "F test" . 

This test compares the relative increase in the residual 
sum of squares by reducing the complexity of the model 
with the relative increase of degrees of freedom, namely, 
the differences between the number of data points and 
the free parameters. Factor F is, thus, computed as 

RS S2 — RSS^ 

n - p3 

where RSS2 and RSS3 are the residual sum of squares 
for the two- and three-phase power-law best fit models, 
P2 — ^ and P3 = 6 are the corresponding number of free 
parameters, and n the total number of points, which, 
in our case, is the effective number of considered mass 
bins (14) and not the total number of variables used in 
the fit (200,000 x 14), given that the additional mul- 
tiplicity has been introduced to map the non Gaussian 
errors. Typically, if the simpler model is valid (which 
means that there is no need to add parameters), then 
F ~ 1, whereas if the complex model represents the dis- 
tribution better, then F > 1. We obtain F = 1.15. The 
probability, P, that the null hypothesis, which is "the im- 
provement found adding additional parameters is solely 
due to chance" (or "the simpler model is good enough" ) 
should be rejected is given by the value of the cumulative 
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Snedecor's F distribution with (p3 — p^^n — ps) degrees 
of freedom. In our case it is P = 0.65, while, generally, a 
value of P > 0.95 is required to state with sufficient sig- 
nificance, that the complex model is better. Under these 
circumstances, we conclude that our system IMF is best 
approximated by a two-phase power-law. 

We compute the uncertainty associated to this result, 
by utilizing a sampling technique. For each mass-bin we 
pick randomly one out of the 200,000 points from the cor- 
responding distribution of ^(Af ) in that bin, and we per- 
form a two-phase power-law fit as described above, ob- 
taining the two measurements for the slope, and the point 
(logM, log^(M)) where the change of the slope occurs. 
After repeating this process for a substantial amount of 
times, we obtain a statistically significant sample for the 
values of the two slopes and the point of change, and 
consequently the relative confidence level of our fit. Fig- 
ure [11] shows the result of our fit to the system IMF of 
LH 95. The best two-phase power-law fit is shown as 
a solid line, while the uncertainty of this fit is demon- 
strated by the grey shaded areas for confidence of 95% 
(dark grey) and 99% (hght grey) respectively. In Table [3] 
we give the system IMF slopes x and their uncertainties 
with the corresponding mass ranges. The la error in the 
determination of the break point is 0.07 dex. 

Another bias that should be considered in the con- 
struction of the system IMF of LH 95 is the possible 
contamination of our selected sample of PMS stars by 
a binary sequence of the LMC field MS population, dis- 
cussed in § 15.11 In order to quantify the effect of this 
bias we select this part of the CMD and remove the cor- 
responding stars. Specifically, we remove all stars bluer 
than the borders defined by the points (F555W— F814W, 
F814W) = (0.9, 22.9), (1.65, 24.3) and (2.2, 26.9) on 
the CMD. The number of these stars changes accord- 
ing to the stochastic Monte Carlo field subtraction; on 
average they count to 82 members, less than 4% of the 
selected PMS population. Then, we reapply the analy- 
sis described in iji^ 15.11 and 15.21 for the construction and 
model- fitting of the system IMF. We find that this new 
IMF docs not differ significantly from the one shown in 
Figure [Til and its shape is not different from that de- 
scribed with the slopes of Table [3] Specifically, this 
IMF has 99%-confidence slopes x ~ 1.16 for —0.5 < 
log(M/Mo) < 0.04, and x ~ 2.03 for log(M/M0) > 
0.04, indistinguishable from the slopes of Table [3] Our 
analysis on the sub-clusters of LH 95 showed that the 
contamination of our sample by a MS binary sequence 
does not alter also the shape of their individual system 
IMFs, with slopes as well indistinguishable from those 
given in Table [H 

The lack of spectroscopy and the use of broad-band 
photometry alone for the mass estimation of massive MS 
stars is found to lead to steeper IMF slopes (|Massev et alj 
|1995[) . and this may affect the overall intermediate- and 
high-mass slope of the IMF of LH 95. H o wever , con- 
sidering that in their analvsis iMassev et al.l (|X995) show 
that this effect is significant for stars more massive than 
15 Mq, the only mass-bin in our IMF that can be bi- 
ased is the most massive one, comprising stars with 
11 < M/Mq < 20. In addition, the IMF derived from 
photometry alone for high-mass stars in star-forming re- 
gions of the Magellanic Clouds is found to be also quite 
shallow and not steep, and within the observed IMF slope 
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Fig. 11.— The system IMF of the association LH 95. The 
best two-phase power-law fit derived in § 15.21 is drawn with a solid 
line, while the shaded areas represent the 95% and 99% confidence 
uncertainties in the slope determination and the break point. Units 
of the IMF are logarithmic number of stars per solar mass per pc^ . 
The corresponding IMF slopes and their uncertainties are given in 
Table [3] 



TABLE 3 

The slope of the system IMF of the 

ASSOCIATION LH 95 



mass range 


X 


0'95% 


"^99% 


-0.5 < log(M/M0) < 0.04 


1.05 


-0.20 
+0.15 


-0.29 
+0.20 


log(A^/M0) > 0.04 


2.05 


-0.28 
+0.39 


-0.32 
+0.53 



Note. — A Salpeter (1955) IMF would have a 
slope of X = 1.35. 



variations in such regions (see discussion in §6.2|1 . Un- 
der these circumstances, the corresponding uncertainty 
in our derived system IMF slope due to the use of pho- 
tometry alone should not be expected to be larger than 
those given in Table 3. 

6. DISCUSSION 

The IMF constructed in the previous section is the 
most complete extragalactic system IMF down to the 
sub-solar regime derived so far. Specifically, our correc- 
tion for incompleteness in the derivation of ^(M) actually 
provides a reliable IMF down to M ~ 0.43Mq, as dis- 
cussed in i i5.1l Probably, the most interesting result con- 
cerning this IMF is the change in its slope, which we ver- 
ified statistically for stars in the the subsolar mass range. 
The IMF slope, x, decreases for these stars by one unit. 
An analogous behavior is generally found in related stud- 
ies in the Galaxy. Comprehensive studies which summa- 
riz e the average properties of the ga l actic IMF a re those 
of 'Scald to86l ). revisited in iScalol (| 19981 ) and iKroupal 
(2001, 2003)- The latter investigations are rather com- 
plete, since they cover stellar mass distributions in a wide 
mass range, from brown dwarfs with M ^ 0.0 1 M(7) t o 
the most massive stars. According to iKroupal (|2002f) . 
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the IMF slope x (Eq. [8]) changes from x — —0.7 in the 
sub-stellar mass range, to a; = 0.3 for masses between 
0.08 Mq and 0.5 Mq, x = 1.3 for 0.5 Mq <M <l Mq, 
and X = 1.7 for stars of higher masses. This IMF is 
generally characterized as the Galactic average, in the 
sense that it is reasonably valid for different regions of 
the Galaxy. A comparison between the system IMF we 
construct here for LH 95 and that o flScaloTdiggl . as well 
as the universal IMF o f Kroupa (|200lD and the average 
Galactic IMF of Krovipa (2002), is shown in Figure [T^ 
We scaled all other IMFs to match our mass distribution 
at log(M/M0) = —0.3. In Figure [12] an overall agree- 
ment between the various IMFs is quite evident. 

The mass limit of the breaking point in our system 
IMF 1 Mq) coincides with the "standard Galactic 
field IMF" , wh ich has been r eported to have a change of 
slope at IMq ()Kroupal l 2002[ ) in ag reement also with the 
knee of the IMF derived bv lScalol (|1998f) . In addition, 
there is a second knee in the Kroupa IMF, which oc- 
curs at the lower mass-limit of 0.5 Mq. Since we cannot 
constrain the system IMF below O.43Af0 we do not iden- 
tify this second knee. The slope of our IMF for stars with 
M ^ 1 Mq falls within the range of slopes found by both 
Scalo and Kroupa. This is in line with the recent finding 
that the system IMF is essentially equal to the stellar 
IMF for intermediate and massive stars (Wcidncr et al. 
l2008l ). Our system IMF, however, is somewhat steeper, 
as an indication of a bottom- heavy IMF, but the use of 
photometry alone for the estimation of the masses for 
the brightest stars (M ^ 15 Mq) may bias the slope 
to steeper values (see iJ5.2p . As far as the slope of our 
IMF in the sub-solar regime is concerned, it is found 
slightly more shallow than that of these authors. Specif- 
ically, our IMF is by Aa; — 0.1 - 0.3 more shallow than 
the stellar IM F (as seen from the theoretical a-plots of 
lKroupal200]l ). making the resulting stellar IMF for LH 95 
indistinguisha ble to the Galactic-field stellar IMF (e.g. 
lKroupall2002[ ). It should be noted that both Scalo and 
Kroupa IMFs extend down to masses as low as the hy- 
drogen burning limit (M ~ 0.08 Mq), well below our 
detection limit. In general, one can conclude that the 
system IMF we derive for LH 95 is in agreement with 
the average galactic IMF in the entire mass range we 
studied, from 0.43 Mq to - 20 Mq. 

The high-mass IMF of LH 95 was constructed earlier 
from direct CCD imaging in BVR w ith the 1-m tele- 
scope at Siding Spring Observatory by iGouliermis et all 
()2002|) . These authors provide the field-subtracted 
completeness-corrected IMF of the system for stars with 
masses down to ~ 2 Mq , and they find that this IMF can 
be fitted with a single power-law with slope x = 1.2 ±0.3 
for stars with masses up to about 17 Mq. This result 
indicates a more shallow slope x than what we found for 
intermediate and high mass stars (x = 2.05). However, 
in our analysis we included masses as low as IMq, ex- 
tendin g the mass range investigated in IGouliermis et ahl 
(I2OOI) . Indeed, if we exclude the mass bins from 1 to 
2 Mq (0 < log (A//Mq) < 0.3) from our analysis we de- 
rive a more shallow value for x for the intermediate- and 
high-mass stars. The effect of a steeper IMF for the mass 
range between 1 and 2 Mq is evident in Figs. [TUl and [TTl 
After excluding stars of M < 2 Mq from the considered 
stellar sample and performing a single power-law fit on 
the resulting £,{M), we find an IMF slope of a; ~ 1.55, 



comparable to that of IGouliermis et ahl (|2002f ). 

Concerning our derived IMF for the central region of 
LH 95 it is important to assess the accuracy of our find- 
ings, in relation with the possible sources of uncertain- 
ties. One of the most important error sources normally 
comes from the incomplete stellar samples. However, in 
our case we consider the completeness as it has been 
derived from artificial star tests in the central part of 
LH 95, where the crowding is higher (Figure HJa), and we 
performed rigorous corrections for incompleteness taking 
into account variations of completeness for all stellar po- 
sitions in the CMD. As far as the very- low mass stars is 
concerned, we do not consider those with M ^ 0.35 Mq, 
where our completeness correction is limited. 

6.1. Stellar Binarity 

An important source of uncertainty in the construc- 
tion of the IMF is also the presence of unresolved bina- 
ries. Unresolved binarity, can modify the shape of the 
observed mass function so that it appears more shal- 
low than what it actually is. For this reason, we stress 
that the IMF reported in Table [3] must be considered as 
the system IMF of our region, and not a stellar IMF as 
which one refers when the effect due to binarity is re- 
moved. Clearly it is not possible to correct this effect 
without the knowledge of the actual binary fraction of 
the studied population, but one can es timate the overall 
effect on the final mass distribution. ISagar fc Richtlerl 
()1991| ) applied Monte Carlo simulations on artificial stel- 
lar populations, in order to investigate the change in the 
mass function slope Aa; due to binarity for stars in the 
mass range 2 < M/Mq < 14. These authors found 
that the result depends on both the binary fraction / 
and the original slope xq of the mass function. Specifi- 
cally, for an orig inal slope x q = 1 .5 and a binary fraction 
/ = 0.5 ISagar fc Rich tleil ()1991f ) found a slope change 
Ax ~ 0.34, while f or xq = 0.5 and / = 0.5, Ax ~ 0.5. 
On the other hand, iKroupa et al.l (|199lD in their investi- 
gation on the low-mass IMF in the solar neighborhood, 
showed that binaries can make a major effect in mass 
function determinations below 1 Mq. Consequently, the 
expected changes in the slope due to unresolved binarity 
may introduce uncertainties in the derived slopes given 
in Table [31 but not large enough to justify the intrinsic 
change of slope that we observe in our IMF for LH 95 
at 1 Mq. In any case, as discussed bv iKroupal ([2(302), 
most of the studied mass functions in the literature are 
not corrected for binarity, mostly due to the high uncer- 
tai nty in the binary fra ction. The Galactic IMF derived 
by (|Kroupall200lLl200^ . however, has been corrected for 
unresolved multiple systems. A more recent study of the 
binary effect on the massive-star IMF bv fWeidner et ahl 
()2008l) shows that a power-law IMF is not significantly 
affected even by large numbers of unresolved binaries or 
higher order multiples. 

In order to analyze the effect of unresolved binarity in 
the measured slope of the system IMF in the sub-solar 
regime, we a pply a Monte-Car l o tech nique adopted from 
the work of ISagar fc Richtled (|1991h and iKroupa et al.l 
(|1991lll993D . For this method we assume a pre-defined 
set of functions for the IMF, and in order to be consis- 
tent with those we find for the sys t em IM F, as well as 
with the Galactic IMF of iKroupal (|2002D . we consider 
three-phase power laws with break points at 0.3 Mq and 
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Fig. 12. — Comparison be tween the deri ved LH 95 system IMF and the average IMF of lScalol 111998') (left). IKroupal II2001I') (middle), as 
well as the universal IMF of lKroupal II2002I ') (right). 



1 Mq. We limit our mass samples in the stellar regime 
with M > 0.08 Mq. We let the slopes xi, X2, X3 within 
the three selected mass ranges to change respectively in 
the range -1 < xi < 1.5 for 0.08 < M/Mq < 0.3, 
< 2:2 < 2 for 0.3 < M/Mq < 1 and 1 < ^3 < 2 for 
M/Mq > 1. For each of these mass distributions, we 
generate a sample of 50,000 stars positioned in the CMD 
along our 4 Myr PMS isochrone (corresponding to the 
age of the system). We assume binary fraction, /, be- 
tween and 1, and we randomly pair couples of stars 
from the simulated distributions. We then compute the 
new position in the CMD of these binary systems and 
assign new masses as if they were single stars. The out- 
coming mass distributions are fitted with a 3-phase power 
law imposing the same position of the break points. 

We isolate the models for which the derived system 
IMF in the range 0.3 < M/Mq < 1 is compatible with 
the one measured for LH 95 within the confidence in- 
terval of the latter. Naturally there are multiple so- 
lutions, given the lack of knowledge about the actual 
IMF of LH 95, its shape at lower masses, and the bina- 
rity fraction appropriate for this system. However, we 
can constrain the uncertainties due to unresolved bina- 
rity in our measured slope, X2, for the subsolar regime 
0.3 < M/Mq < 1. The IMF slope within this mass range 
is mainly affected by the value of /. If we assume xi to be 
that of the Galactic-field IMF, the slope X2 in this mass 
range changes by Ax2 — 0.45 for / = 1. This change 
decreases to Ax2 — 0.35 for f = 0.6, and Ax2 — 0.25 for 
/ = 0.4. 

A change in the original IMF for masses larger than So- 
lar does not affect our results by more than Ax^ = 0.1, 
whereas changing the slope at very low masses (0.08 < 
M/Mq < 0.3) affects significantly the slope X2, diminish- 
ing it when xi becomes much steeper or shallower than 
the value of lKroupal ((2002) IMF. Even in this case, the 
effect is always less than Ax^ = 0.2. Therefore, we con- 
clude that for a reasonable binary fraction of / ~ 0.5, 
unresolved binarity can bias our IMF slope, x, in the 
subsolar regime on average by no more than 0.3 units. 
As a consequence unresolved binarity does not affect sig- 
nificantly the measured IMF for LH 95 within the es- 
timated errors, and thus the corrected stellar IMF for 
LH 95 remains indistinguishable from the Galactic IMF. 



6.2. Metallicity effects 

Recent findings suggest that variations in the metal- 
licity of a star-forming cloud might cause variations in 
the IMF slope for stars with masses M < 0.7 Mq, in the 
sense that Galactic regions with higher \Fe/H] appear 
to produce more low-mass stars ()Kroupal |2002[ ). This 
is further support ed by the measur ed IMF in Galactic 
open cluster s (Barrado v Navascues et aLll2001[ ). globu- 
lar clusters (jPiotto fc Zoccalil I1999D and old and rela - 
tively rnetal-p oor thick disk stars (|Revle fc Robinll200l() . 
iKroupal (|2002l ) proposes a systematic metallicity depen- 
dence of the IMF exponent of the form x w 0.3-1- 
Ax[Fe/ H], where Ax « 0.5. If we apply this de- 
pendence to the LMC, for which [Fe/H] ~ —0.4, we 
find that the slope of the average IMF in the LMC for 
masses M < 0.7 Mq should be fiattened to the value 
X « 0.1. The flattening of the IMF of the LMC general 
field for stars with AI < 0.7 Mq has been discussed by 
iGouliermis et al.1 (|2006al ). who found only indications of 
such a change in the IMF slope, because of observational 
limitations. Our data, on the other hand, being deeper, 
show a constant IMF slope in the low-mass regime down 
to the detection limit of M ~ 0.43 M© with no indication 
of fiattening for M < 0.7 Mq. Consequently, we cannot 
verify any metallicity dependence of the sub-so lar IMF 
in the LMC, at least based on the formula of IKroupal 
(|2002f) . Considering that this result is based on a single 
stellar system, more studies of low-mass populations in 
the Magellanic Clouds may provide a better constraint on 
any systematic metallicity-dependence of the low-mass 
IMF. 

As far as the high-mass regime is concerned, while the 
IMF of LH 95 seems to be steeper than that derived 
from ground-based observations by IGouliermis et ahl 
(|2002l) . our analysis (fj6]) shows that if we consider 
only stars more massive than 3 Mq the IMF of 
LH 95 has a Salpeter-li ke slope, comparable to that of 
IGouliermis et al.l ()2002D . While, as discussed in ^5.21 
the use of photometry alone leads to steeper IMFs, 
there are several photometric studies of associations in 
the LMC, which derive a high-mass IMF that is not 
steeper, but more shallow than the typical Salpeter 
IMF (e.g. iHill et al.lll99l l()^[T99^ iDolohin fc H"u^ 
I1998D . Taking into account all available studies for the 
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Fig. 13. — Measured mass function for the 3 subclusters of LH 95, 
as they are defined in Table[2]and shown in Figure|4}3. The method 
used for the determination of the best fitted two-phase power-law 
is analogous to that applied for the entire central region of LH 95, 
except for the fact that here we constrain the break point to be the 
same derived for the global IMF, at log(M) = 0.04. The overall 
§(M) of LH 95 is overlayed as a dashed line. 



intermediate- and high-mass stars with M > 7 M 



both s pe ctroscopic and ph otometric (e .g. 'Massei^ 



Hi 



, et aLl 

1989allbl: iHill et al.iri 994: Massev et all Il995: Hil l et all 
1995t Massev fc Huntenil998f ) one can see that the corre- 
sponding IMF slopes in star-forming regions of the LMC 
are found clustered around the value x « 1.5, similar 
to the Galactic average IMF. The system IMF we de- 
rive here for LH 95 for stars within the same mass limits 
does not differ significantly from this value. Extensive 
studies have also shown that the high-mass IMF is inde- 
pendent of metallicit y in the Galax y and the Magellanic 
Clouds (see review bv lMassevl[2003l ). and our results also 
suggest a weak metallicity dependence of the high-mass 
IMF, providing further support to its universality, in the 
sense that the average IMF of a galaxy seems to be inde- 
pendent from environmental factors, and that any indi- 
vidual stellar system is not expected to share the exact 
IMF Indeed, the stan dard gal actic IMF, as presented by, 
e.g.. [ScaIol (pM ) and lKroupal ((200211 . is the average over 
a set of systems, each of them showing slightly different 
mass distributions. 

6.3. Variability of the IMF within LH 95 
The clustering behavior of the PMS stars in LH 95 



is discussed in § 13. 2| where three high concentrations 
of such stars within the main body of the system are 
identified as the sub-clusters of the association LH 95 
(given the names cluster A, B and C respectively). In 
order to study how the IMF may vary from one sub- 
cluster to the other, and how the IMF of each individual 
concentration of PMS compares to the "overall" IMF of 



TABLE 4 

IMF SLOPES AND UNCERTAINTIES FOR THE THREE SUB-CLUSTERS 

OF LH 95 




-0.5 < log(A^/M0) < 0.04 


log{M/MQ) > 0.04 


subcluster A 


0.94;«:« 




subcluster B 




1 ao-0-60 
^•^^+1.71 


subcluster C 


1 

'^■"^+0.33 


1 70-0.29 



Note. — slopes are in terms of x according to Equation \8\ 
uncertainties are 95% confidence intervals. 

the system, we constructed the IMF for each of the three 
sub-clusters. 

We applied the same method used in i i5.1l and l5.21 iso- 
lating each time the regions of the sub-clusters, as shown 
in Figuredb, iterating the Monte Carlo field subtraction, 
computing the distribution of number of stars in mass 
bins and correcting for incompleteness. For the model fit- 
ting, we used the same Levenberg-Marquard technique, 
with the diff'erence that, in this case, we constrained 
the position of the break point between the two power- 
laws to be the same as for the general IMF, namely, at 
log(Af/AfQ) — 0.04. This choice is reasonable consider- 
ing that we are mainly interested in eventual changes in 
the derived slopes of the IMF. 

The constructed IMF for each sub-cluster is shown in 
Figure [m The overall IMF for the whole system is also 
overlayed, after it is normalized in order to match the 
subclusters' best fit on the break point. In each case the 
IMF for the subsolar masses is the same as the general 
IMF, as discussed in ^for the entire system of LH 95. 
Moreover, the derived two-phase models that fit the IMF 
of each sub-cluster and the corresponding uncertainties, 
given in Table [H are compatible to that of the overall 
IMF given in Table El 

The IMF constructed for each sub-cluster in LH 95 sug- 
gests that the mass distribution of the PMS stars in this 
association does not vary much with location. Moreover, 
the spatial distribution of the intermediate- and high- 
mass stars within the whole of the system does not show 
any indication that the system is mass segregated. On 
the other hand if the individual sub-clusters are consid- 
ered, the more massive stars do occupy the central parts 
of each of them, providing clear evidence that mass seg- 
regation in LH 95 probably occurs in rather small clus- 
tering scales. 

Acc ording to recent results by IWeidner fc Kroupal 
(|200(l) on the integrated galaxial stellar IMF, the over- 
all IMF of a composite system should be steeper than 
the one it would have if it was not composite, in the 
intermediate- and high-mass regime. Taking into account 
the fact that LH 95 is resolved into at least three individ- 
ual subclusters fi i3.2p . it may be considered as a scaled- 
down case of the scenario of these authors, and therefore 
as a composite system. It should be noted that also this 
may explain the somewhat steep slope of the IMF we 
derive for LH 95. 

7. CONCLUSIONS 

We present a comprehensive investigation of the stel- 
lar populations of the star-forming association LH 95, lo- 
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cated on the north-western edge of the super-giant sheU 
LMC 4 in the Large Magellanic Cloud, using of deep 
HST/ACS photometry in the filters F555W^ and F814W^ 
(proxies for standard V and / respectively). Our findings 
can be summarized as follows. 

We construct a new set of observational PMS mod- 
els (both isochrones and evolutionar y tracks) by start- 
ing from the original computations bv lSiess "etall (iOOOh 
and performing synthetic photometry on artificial stellar 
spectra. The latter were derived from interpolatio n of at- 
mosp here models fr om the Next Gen (Hau sch ildt et al.l 
IT999h and Kurucz (jl uruc j [1993I ) grids. This approach 
enabled several improvements considering observed data 
for PMS stars. Specifically: 

(i) We computed absolute magnitudes for each point 
of the original theoretical models in 88 different 
photometric bands, including the ACS/WFC sys- 
tem, which is useful for our observations. The avail- 
ability of observational PMS evolutionary models 
in any given photometric system allows us to avoid 
the application of any color corrections to the ob- 
served data, or to approximate observations in the 
ACS / WFC system with any standard photometric 
system. 

(ii) Our com putations, being similar to that of 
iGirardi et al. (2002), do not rely on a single de- 
pendence on effective temperature for the compu- 
tation of colors and magnitudes, but they also take 
into account the dependencies on the stellar surface 
gravities (i.e. in luminosity), and therefore they 
provide a more accurate treatment of the young 
stages of PMS evolution. 

(iii) These observational models ar e extrapolated to 
lower metallicities than those of iSiess et al.l ()2000f ) 
by computing our synthetic photometry for the av- 
erage LMC and SMC metaUicities of Z = 0.008 
and Z = 0.004, respectively. 

We enhanced the catalog of stars detected in the ob- 
served field around th e association LH 95 with ACS 
(|Gouliermis et al.l2007( ) by identifying new point sources, 
which were previously rejected by the photometry mostly 
in crowded regions. Therefore, about 900 additional stars 
are included in our photometric catalog (Table 1), which 
now includes in total 17,245 stars. Our deep ACS pho- 
tometry allowed us to analyze the properties of the abun- 
dant PMS population found in LH 95, which shows an 
average age of 4 Myr. We performed star counts on our 
stellar sample of the area of the association and we se- 
lected the central part of the observed field, where a sta- 
tistically significant concentration of PMS stars is found, 
as the most representative region of LH 95. The Scr iso- 
pleth was selected to define the boundaries of the sys- 
tem. The PMS stars within the association are found to 
be mostly concentrated in three distinctive sub-clusters, 
named cluster A, B and C respective, providing evidence 
of the existence of stellar subgroups within a single asso- 
ciation. 



We estimate the total stellar masses and corresponding 
dynamical timescales for the system itself as well as the 
three subclusters, and we find that the system cannot be 
considered dynamically evolved, and it does not experi- 
ence any disruption. Moreover, the fact that LH 95 is a 
bound cluster with its gas partially removed implies that 
the star formation efficiency is probably higher than the 
typical values for young clusters. 

We estimated the stellar contribution of the general 
background field of the galaxy in the area of the system 
from our ACS observations of a nearby control field. We 
then performed field subtraction from the area of the 
system with the use of Monte Carlo simulations. 

With our evolutionary models we estimated the inter- 
stellar extinction in LH 95 and we computed the masses 
of all the members of the system. Our stellar sam- 
ple reaches masses as low as ~ 0.2 Mq, and we are 
able to correct the incompleteness accurately down to 
~ 0.43 Mq. For the first time, this allows the construc- 
tion of the extragalactic IMF down to the sub-solar mass 
regime. 

We derive the system IMF by counting stars within 
bins of variable width, which increase with mass, and 
correcting for incompleteness based on artificial star ex- 
periments for the stars included in each mass range. We 
derive the actual statistical distribution of the counting 
errors associated to each mass bin via a Monte Carlo 
technique, in order to take into account uncertainties in- 
troduced by the contamination from the field population. 
The statistical F test for model selection allowed us to de- 
termine a two-phase power-law as the most appropriate 
model to fit our IMF. We detect a change in the slope 
of the IMF (the "knee" of the IMF) for masses lower 
than 1 M0. A Levenberg-Marquard algorithm for least 
square fitting was applied to determine the power-law 
exponents and the point where the slope changes. The 
derived system IMF slope is found to be a; = I.OS^q;^^ in 
the subsolar regime, and x = 2.05^g'3g at intermediate 
and high masses, with a 95% of confidence. 

This system IMF, once corrected for unresolved bina- 
rity, is statistically compatible, within the studied mass 
range, with the average Galactic IMF, concerning its 
slope and the mass limit where it changes (1 M0). As 
a consequence, we cannot consider the metallicity differ- 
ence between the LMC {Z = 0.008) and the Galaxy im- 
portant enough for a significant change in the IMF slope. 
Certainly similar data on such stellar systems will show 
if this is a normal behavior of the IMF in this galaxy. We 
find no significant differences in the shape of the overall 
IMF of LH 95 from that of each of the three individ- 
ual PMS sub-clusters of the association. This clearly 
suggests that the IMF of LH 95 is not subject to local 
variability. 
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